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By devising a novel framework, we present a comprehensive theoretical study of solubilities of 
alkali (Li, Na, K, Rb, Cs) and alkaline earth (Be, Ca, Sr, Ba) metals in the he boron-rich Mg-B 
system. The study is based on first-principle calculations of solutes formation energies in MgB2, 
MgB4, MgB7 alloys and subsequent statistical-thermodynamical evaluation of solubilities. The 
advantage of the approach consists in considering all the known phase boundaries in the ternary 
phase diagram. Substitutional Na, Ca, and Li demonstrate the largest solubilities, and Na has the 
highest (0.5-1 % in MgB 7 at T = 650 - 1000 K). All the considered interstitials have negligible 
solubilities. The solubility of Be in MgB7 can not be determined because the corresponding low- 
solubility formation energy is negative indicating the existence of an unknown ternary ground state. 
We have performed a high-throughput search of ground states in binary Mg-B, Mg-A, and B-A 
systems, and we construct the ternary phase diagrams of Mg-B- A alloys based on the stable binary 
phases. Despite its high temperature observations, we find that SrgMg38 is not a low-temperature 
equilibrium structure. We also determine two new possible ground states CaB4 and RbB4, not yet 
observed experimentally. 

PACS numbers: 74.70.Ad, 61.50.Ks, 81.30.Hd 



I. INTRODUCTION 



diminished by a segregation into LiB phase. 



The interest in magnesium diboride emerged after the 
discovery of superconductivity in MgB2 at about T c =39 
Attempts to increase T c by small additions of al- 
kali NaZ£, Rb£±2, Cs&I£) and alkaline earth 
(Beii, Caii±2iH, Ba£) metals to MgB 2 , proved to be un- 
successful. The difficulty was attributed not only to the 
inability of such solutes to decrease T c but also to their 
low solubility and precipitation in secondary phases. Al- 
though a claimed superconductivity at 50K was reported 
for the Mg-B- A (A=Cs, Rb, Ba) system^, attempts to 
reproduce the results have been so far unsuccessful (for 
A=Cs, Rb)i^. The problem can be attributed to solutes' 
segregation in grain boundaries and to thus to low sol- 
ubility in the bulk phase (much lower than reported in 
Ref. i). 

For consistent interpretation of experimental observa- 
tions, theoretical studies of solubility in Mg-B of alkali 
metals and alkaline earth metal solubilities are therefore 
necessary. For instance, in Ref. |~L4|, the semi-empirical 
Miedema approach-^and the Toop's models were used 
to address the the heats of formation of binary alloys Mg- 
A, B-A, and Mg-B (A=Li, Na, and Ca). It was proposed 
that Ca may form stable compounds while Na, Li could 
lead to meta-stable or unstable ternary phases in MgB2. 
In Ref. UtI by calculating first-principle formation ener- 
gies of Li and Na impurities in MgB2, and by neglecting 
the effects of other ground states in the ternary phase 
diagram solubility calculations, it was concluded that Na 
should have very low solubility, whereas the solubility 
of Li should be comparatively higher although modestly 



The present paper is orthogonal to previous studied. 
We develop a comprehensive theoretical framework to 
determine the solubilities of alkali metals (Li, Na, K, Rb, 
Cs) and alkaline earths (Be, Ca, Sr, Ba) in the boron-rich 
Mg-B system. The study consists of first-principle cal- 
culations of solutes formation energies in MgB2, MgB<4, 
MgBy alloys and subsequent statistical-thermodynamical 
evaluation of solubilities with respect to all known equi- 
librium states of the Mg-B-A system. The results help 
outlining future directions in experimental searches. 



The paper is organized as following. In Section [U we 
describe the adopted solubility mechanisms in Mg-B. In 
Section lLnl we introduce the relevant impurity formation 
energies in terms of supercell energy calculations and the 
appropriate ground state(s). In Section llVl the approx- 
imation for the free energy of Mg-B-A solid solution is 
formulated. In Section |Vl we present an approach for 
solubility calculation considering all the ternary ground 
states. A simple analytical low-solubility approximation 
is devised. Section IVT1 is devoted to the high-throughput 
ah initio search for ground states in binary Mg-B, Mg- 
A, and B-A systems, to the ternary phase diagrams of 
Mg-B-A systems, and to the impurity formation energies 
determined through the phase boundaries of the systems. 
The numerical values of solubilities are presented in Sec- 
tion lVIII Scction lVIIII summarizes the results, draws con- 
clusions, and comments on strategies for future research 
in this area. 
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II. MG-B-A SOLID SOLUTION 

For the description of solubility of alkali and al- 
kaline earth metal elements A (A=Li, Be, Na, Mg, 
K, Ca, Rb, Sr, Cs, Ba) in the Mg-B system, we 
consider the disordered solid solutions of A-atoms 
as interstitial and magncsium-substitutional impurities 
in the experimentally reported compound a 18 ' 19 MgB2, 
MgB4, and MgBy. We do not consider boron sub- 
stitutions by A because they has not been observed 
experimentally^^^^ii 2 .^. The disordered 
solid solution of A inside Mg-B is labeled as "(1)" 
throughout the paper. 

While the magnesium-substitutional positions are 
determined^, the "most accommodating" interstitial lo- 
cations have to be found. The task is implemented with 
the following exhaustive search performed through our 
software AFLO W 20 ' 21 . Let us consider a quadruplet of 
no-coplanar atoms, where the first atom belongs to the 
unit cell and the others are closer than the maximum 
diagonal of the unit cell to the first atom. A cage is 
defined when the spherical region of space touching all 
the four atoms of the quadruplet docs not contain fur- 
ther atoms inside. An interstitial position is found if the 
cage has its center inside the unit cell. By considering all 
the possible combinations, the symmetrically incquiva- 
lent interstitials can be identified through the calculation 
of their site symmetry (with the factor group of the unit 
cell). Note that in unit cells with complex arrangements, 
many of the interstitials positions can be extremely close. 
Thus, an interstitial atom located in any of those close 
positions would deform the nearby local atomic environ- 
ment and relax to the same final location. Hence, the 
number of symmetrically inequivalent cages can be fur- 
ther reduced by considering the whole set of positions 
that would agglomerate upon insertion of an interstitial 
atom, as a single interstitial position. The results of the 
search arc presented in Table [U The table demonstrates 
that the higher boron contents the larger number of cages 
and the larger radius of the bigger cage. 



III. FORMATION ENERGIES DEFINITIONS 

"Raw" formation energies 

Let us define the so-called "raw" formation energies 
Hi and ~p s (composition unpreserving-2 2 -) as the changes 
of the energies of the solvent upon introduction of one 
solute atom A in the i-th type interstitial or s-th type 
substitutional positions. In first-principles calculations, 
the solvent can be replaced by a large supercell ( "sc"), 



TABLE I: Centers and sizes of the symmetrically inequivalent 
interstitial cages for the MgB2, MgB4, and MgB7 unit cells. 
Coordinates are presented as fractions of the standard unit 
cell basis vectors (a, b, c)^. The radius is the distance from 
the cage center to the nearest Mg or B atom. Vi indicates the 
numbers of the symmetrically equivalent cages within the unit 
cell. The total numbers of distinct geometrical cages before 
agglomeration upon insertion of an interstitial is reported in 
brackets. 



Structure 


i 


Coordinates 


Cage 










(fract. of a 


,b,c) 


radius (A) 




Mj 


5B 2 


1 


0.341, 


0.681, 


0.000 


1.7619 


2 (6) 






2 


0.006, 


0.000, 


0.500 


1.7616 


1 (6) 






3 


0.500, 


0.500, 


0.122 


1.6015 


6 (6) 


M( 


?B 4 
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0.106, 


0.669, 


0.412 


1.9145 


4 (52) 


(Mg 


iBifj) 
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0.375, 


0.250, 


0.337 


1.7836 


4 (4) 






3 


0.145, 


0.632, 


0.813 


1.7821 
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0.144, 


0.750, 


0.925 


1.7252 


4 (24) 
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0.225, 


0.750, 


0.719 


1.6322 


4(4) 
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0.086, 


0.015, 


0.527 


1.6128 


8(8) 
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0.101, 


0.543, 


0.261 


1.5322 


8(8) 


Mj 


?B 7 
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0.000, 


0.049, 


0.618 


2.0050 


8 (152) 


(Mg 8 B 56 ) 
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0.189 


0.250 


0.821 


1.8937 


8 (24) 
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0.250, 


0.217, 


0.750 


1.7995 


8(8) 
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0.243, 


0.250, 


0.505 


1.7888 


8 (24) 
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0.197, 


0.250, 


0.469 


1.7802 


8 (24) 
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0.000, 


0.250, 


0.150 


1.7002 


4 (24) 






7 


0.012, 


0.000, 


0.000 


1.6995 


8 (64) 






8 


0.250, 


0.062, 


0.250 


1.6068 


8 (40) 






9 


0.000, 


0.171, 


0.788 


1.4991 


8(8) 






10 


0.131, 


0.186, 


0.119 


1.4443 


16(16) 



so that 

Hi = 4 ( A ) = £sc[A (l) Mg, isc B„ l3c ] - E BC [H Sns B m J, 

7 i sS 4 1) (a)-^ 1) (^)= (1) 

SsctA^Mg^.xBmJ - E sc [n &ns B m J, 

where n BC and m BC are the numbers of Mg and B atoms in 
the supercell; ^'(Mg), ^ (a), and /is (a) are the zero- 
temperature chemical potentials of the Mg-atom, and the 
A-atoms in the i-ih interstitial and s-th substitutional 
positions, respectively 2 ^. 

The conventional unit cells of MgB2, MgB/i, and MgBy 
are used to construct the appropriate supercells. The 
supercells are chosen to have solute-solute A- A distance 
at least ^2.5 times the nearest neighbor solvent bonds in 
order to diminish the contribution of the A- A interactions 
to the calculated energies. We use 2x2x2 unit cells for 
MgB2, while there is no need to create supercells of MgB4 
and MgBy since their unit cells arc already large enough. 
Note that boron-boron (B-B) bond is the shortest among 
the Mg-Mg, B-B, and Mg-B bonds. Table IITI summarizes 
the parameters of the supercells and Table IIIII lists the 
supercell atomic compositions. 
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TABLE II: Numbers of atoms n,m (Mg n B m ) in the con- 
ventional unit cells ( "uc")^ and in the constructed supercells 
( "sc") for MgB2, MgB.4, and MgBy. v a indicates the number 
of symmetrically equivalent s-type Mg atom positions within 
the corresponding unit cell: ^ s v s — n U c- d^ and are 
the distances between the nearest solutes and nearest solvent 
atoms (boron-boron in our case) in the supercell. 



Compound 


nnc 






rise 


m S c 


d B B (A) 


dkk/d BB 


MgB 2 


1 


2 


1 


8 


16 


1.7811 


3.46 


MgB 4 


4 


16 


4 


4 


16 


1.7004 


2.60 


MgB 7 


8 


56 


4,4 


8 


56 


1.7372 


3.44 



TABLE III: Interstitial and substitutional supercell atomic 
compositions (a;| c ,2:| c ,2:^g) as functions of the nsc (for Mg) 
and m S c (for B) numbers of atoms in the super cell. 





Interstitials 


Substitutional 


~.sc 

x k 

J,g 
~sc 


[1 + rise + m S c]~ L 
m sc [l + rise + msc] -1 
n S c[l + rise + rase] -1 


[rise + msc]~ L 
m sc [n S c + m S c] _1 
n S c[n S c + msc]" 1 



"True" formation energies 

The "true" formation energies (composition 
preserving^) are defined from the "raw" formation 
energies (composition unpreserving^) as 



E { ^ = fM - fi^ 



pure 



pure 



pure 

Ma . 



(2) 



where ^^ re and [J^^ 6 arc the chemical potentials (en- 
ergies per atom) of A- and Mg-atoms in pure A- and 
Mg-solids at zero temperature, respectively. 



Low-solubility formation energies 

In Ref. [24|, a quantity called "low-solubility formation 
energy" , was shown to determine the dilute solubility in 
binary alloys through the temperature exponential fac- 
tor. The quantity can be generalized to the case of Mg- 
B-A alloy as (compare with Eq. (29) in Ref. [13): 



E: ol (A) ee AE^Jxf 
with a = (i), (s) and 

E S c[A (i) Mg n B ms< 



(3) 



AE. 



AE. 



(i) 

sc/at l+nsc+™sc 
£sc[A (fl) Mg B„ 



sc/at 



"sc+msc 



pmix/ sc sc sc\ 
-^at V X A i X B > x Mg/>' 

_ cmix/ sc sc „,sc\ 
■^at V X A ' X B ' x Mg^' 



where E^ x (xf c , x% c , x^) is the energy of a the three 
phase mixture at the same composition of the supercell 
(see Table [m|: 

pmix/ sc sc sc\ _ Y -^uc [Mg„ ac B„ 1|ic ] ( 2 ) (3) 



l Mg> - A l - +m T^i^t +^3^ at 

(5) 

The coefficients Xk (£;=1,2,3) are determined from the 
following linear system of equations 






X 2 xf - 


f X 3 xf ] 




?n uc 


+ X 2 x^ ] - 


f X 3 x ( ^ 


(6) 


^UC 


+ * 2 4g - 







(4) 



where ^at'^A^B (fc=2,3) are the energies (per 

atom) and stoichiometric compositions of the two other 
ground states that, together with Mg niic B muc , form the 
convex-hull triangle containing the point (xf c , x§ c , xfjg) 
in the ternary phase diagram. For example, MgB2, 
Mg, and LiB 3 are the three ground states surrounding 
Mg7Bi 6 Li (a supercell of MgB 2 with a substitutional Li 
atom) as shown in Fig. [5] The coefficient Xk (A:=l,2,3) 
represents the fraction of the k-th phase in the mixture. 
The quantities AE^^^ defined in Eq. |[4j are the su- 
percell formation energies (per atom) determined with 
respect to the mixture of ground states^. E£ o1 (A) and 
AE™ C , at can not be negative. If they do, it indicates that 
the list of ground states is incomplete, and a better phase 
diagram should be established (the missed ground state 
can be supercell itself)^. 

In analogy of with binary alloys^, the ternary alloy 
"low-solubility formation energy" E£ 0l (A) is shown to 
determine the solubility of A in the low-solubility limit 
(Section E| . 



IV. FREE ENERGY OF MG-B-,4 SOLID 
SOLUTION 

The Gibbs free energy per unit cell ( "uc") of the Mg- 
B-A solid solution is determined within the mean-field 
approximation: 

Guc[{ci, c s }, T] = £ uc [M gnuc B muc ] + Aft + A 5s , 

Aft = J2i Viici^i + ^bE[c % Incj + (1 - Ci)ln(l - a)]}, 

Ag s = McsK + k B T[c s lnc s + (1 - c s ) ln(l - c s )]}, 

(7) 

where fcs is a Boltzmann constant, T is the temperature, 
and E uc [w.g nuc B muc ] represents the energy (enthalpy) of ini- 
tial Mg n B TOuc unit cell without A-solutes; the summations 
over i and s arc over all the inequivalcnt types of inter- 
stitial and substitutional positions in unit cell; Ci and c s 
are equal to the site-concentrations of A-atoms at each 
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interstitial (i-) and substitutional (s-) type, respectively. 
Thus, we assume that the concentrations in the equiv- 
alent positions are equal in the disordered state. The 
"raw" formation energies \ii and JL S are introduced in Sec. 
Hill Although, the mean-field approximation neglects 
correlations, it should work well when the deviation from 
stoichiometry is small (see Sec. 19 in Ref. l25l ). In ad- 
dition, we neglect solute-solute interactions that might 
be important especially for high solute concentrations. 
In conclusion, our model is similar to Wagner-Schottky 
model of a system of non-interacting particles^. 

For given concentrations {ci,c s }, the total number of 
atoms per unit cell, N**, the total concentrations of A- 
interstitials, xj^]), A-substitutional, x^ s) , and the total 



concentration of A-atoms, , RFC determined as 



7V at 

1 "uc 



at 
uc ! 



(8) 



'am 



C A ( 3 )- 



At given temperature and concentration , the 

Gibbs free energy per atom G^.' \ x l i T\ is determined 
by minimizating Eq. ([7]) with respect to {c,} and {c s }: 



G£i4\T} = mm G^[{ Cl ,c s },T] 



(i) ' 



where 



G« [{ Cl , c s } , T] = G & [{c, , c s } , T] /AT, 



(9) 



(10) 



The solution {ci,c s } defines the equilibrium distribu- 
tion of interstitial and substitutional A-solutes in the 

M gnuc B m uc Solvent. 

In the case of small concentrations of interstitials, the 
minimization can be done with the Lagrange multiplier 
method, obtaining: 



1 + cxp - — 



k B T 



1 + cxp 



k B T 



(11) 

where the Lagrange multiplier /j, is determined from the 
following equation (derived from Eq. ©): 



x^\n uc + m uc ) = ^ ^Cj(l - ) + 



can differ from stoichiometry through solution because 
of the entropic promotion (Sec. IVIII and Ref. l24l ). At a 
given temperature, the solubility of A-atoms in a com- 
pound is defined as the maximum homogeneously achiev- 
able concentration of ^4-atoms, without the formation of 
a new phase. 

To calculate the solubility- 2 ^, we consider the Gibbs free 
energy G™£ x of the mixture of three phases with a given 
general ("gen") composition x gen = (xf en , x\ en , x^ 1 ) . 
The first phase is the substitutional and/or interstitial 
solid solution with atomic concentration x^ of element 
A in Mg„ uc B muc . At zero temperature, the two other 
phases and Mg„ uc B TOUC form a triangle containing the 
point x gen in the ternary phase diagram. At a finite 
temperature, the free energy G™^ x is the generalization 
ofEq. ©: 



X 2 [x£\T\E$+X a [x)t\T\E$>, 



(1) 



A3) 



(13) 



where G$[x%\T\ is given by Eq. ©, and the second 
and third phases are assumed to be stoichiometric, so 
that their free energies can be approximated by their 
ground state energies (enthalpies) E^ t and E^ t , respec- 
tively. The approximation does not affect much the re- 
sults, as we are interested in the small solubility regime of 
A. Similarly to Eq. ([6]), the fractions Xk are determined 
by solving the system: 



„s en _ 



„s en _ 



= X x x 
Xix- 



(i) 

A 

(1) 



gen _ 



= XiX- 



(1) 
Mg 



X 2 x- 



(2) 
A 

(2) 
B 

(2) 
Mg ■ 



X 3 x 



(3) 
A 
(3) 
B 

(3) 
Mg i 



(14) 



where the second "(2)" and third "(3)" phases are at 
stoichiometry. In addition, we have 



_ "uc 



"HC+niuc 
_ 1 



1 



„(!) 



'Mg 



x {l) \x {1) T 



(15) 



where 



A ' 



T] is the equilibrium concentration of 
A- interstitials in the first phase given by Eqs. (JSJ9J) at 
for chosen total concentration x^ . The minimization of 



G™t x with respect to x^' gives the solubility x\ ' (T) m 



(12) the first phase "(1)": 



V. SOLUBILITY 

According to Nernst's theorem, either a single com- 
pound or a phase separation of compounds at correct sto- 
ichiometry can be present at equilibrium at zero temper- 
ature. At finite temperatures, the composition of phases 



G£f[4 1J (T),T]=minG 



T] 



(i) 

A 



(16) 



This procedure, equivalent to the common-tangent 
method's the generalization of the approach developed 
in Ref. [2J to the case of ternary alloys. The method is 
somehow different from that developed in Ref. 12a. Al- 
though, the regular solution model used in Ref. l28| and 
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the presented ideal solution model coincide in the low 
solubility regime, the consideration of only two ground 
states ( "gs") by the authors of Ref. HH, differs from our 
approach requiring the knowledge of the whole ternary 
stability. This is because the disordered Mg-B-A solu- 
tion does not generally belong to the Mg-B<->gs( 2 ) or 
Mg-B<-> gs^ lines in the ternary phase diagram. Thus, 
we have to consider the phase mixture of Mg„ uc B„ luc , 
gs( 2 \ and , gs^ to guarantee accurate estimation of the 
solubility. 



Low-solubility approximation 

In order to get the analytical expression for equilib- 
rium solubilities from Eq. (|16p . further approximations 
are required: (a) the equilibrium concentration of solute 
A is small and (b) only substitutional or interstitial posi- 
tions of one type a are occupied (this approximation will 
be eliminated at the end of section). Thus, from Eqs. 
(113ll6j) ) we obtain: 



dx 



(i) 



dx 



(i) 



T] 



dx 



(i) 



0, 



where 



KB® 



A 2 ) 

J SLt 



(17) 



(18) 



and the fractions X k = Xk[x^ , x gen ] come from Eqs. 
(|14I15|) . Note that the circa («) in Eq. (fT7|) corresponds 
to approximation (a) applied to Eqs. ()7ll0p : 



dX x (1) 



dx 



(i) 



dx 



(!) ^at 



Equations (I14I15[) lead to: 
dX k 



dx^ 

and through Eq. (fT8| : 



(i) 



dx 



(i) 



dX k 

d xr 



dxr 



(19) 



(20) 



(21) 



By using the expression of G' 1 ' from Eq. ([7]), we have 



dG 



(i) 



dx 



(i) 



Ha + k B Th\ 



1 



(22) 



where the circa («) corresponds to approximation (a)^. 

Equation (jTTJ) can be solved with respect to c a with 
the help of Eqs. (|21l22p as 



1 + exp 



£ S Q ol(^) 



exp 



k B T 



(23) 



'=B T « i! s i(^) 



where 



E% ol (A)=» a - 



dxr 



„s en 



E. 



(i) 



gen 



-IE, 



a£/ at gen 

dx? n 



e£\ x z™) 



^mix 



(xS en ) 



(24) 



Both Equations (|24p and © define the same quantity: 
the low-solubility formation energy, E^ o1 (A), determin- 
ing the low-solubility of A-solutcs in a type positions. In 
addition, by using Eqs. ([8|)- ([23|l ^ the total equilibrium 
concentration of solute A in the phase (1) becomes 



(25) 



where c Q is determined by Eq. (|23[). 

Approximation (b) can be relaxed if the various types 
of substitutional and interstitial positions are occupied 
independently. This is expected to be true in the low 
solute concentration limit. Thus, expression (|25[) can 
be integrated out through the various types of positions, 
leading to 

(1) _ E l v jc± + E s ^£f 



(26) 



where Ci and c s are obtained from a set of Eq. (|23[) , by 
using the corresponding E^^A) and E^^A) for each 
type of defects. 



VI. FIRST-PRINCIPLES CALCULATIONS 

The first-principles calculations of energies are per- 
formed by using our high-throughput quantum calcu- 
lations framework AFLO W 20 ' 30 i 31 ' 32 and the software 
VASRH. We use projector augmented waves (PAW) 
pscudopotcntials 3 ^ and exchange-correlation functionals 
as parameterized by Perdew and Wang 3 ^ for the gener- 
alized gradient approximation (GGA). Simulations are 
carried out without spin polarization (not required for 
the elements under investigation), at zero temperature, 
and without zero-point motion. All structures are fully 
relaxed (shape and volume of the cell and internal po- 
sitions of the atoms). The effect of lattice vibrations 
is omitted. Numerical convergence to within about 1 
meV/atom is ensured by enforcing a high energy cut-off 
(414 eV) and dense 4,500 k-point meshes. 



Ground states determination 

The calculation of solubility of A in Mg-B compounds 
requires the knowledge of the relevant ground states in 
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the ternary Mg-B- A system. The systems under investi- 
gations have not been well characterized experimentally 
or theoretically, and only the binary Mg-B, B-A, and Mg- 
A systems have been studied. Hence, we performed ad- 
ditional high-throughput searches to determine if further 
ground states exist in the three binary system o 20 ' 30 ' 31 ' 32 . 
Based on the knowledge of the binary systems, we buillj 3 ^ 
the ternary ground state phase diagrams for Mg-B-A, 
with the expectation that no missed ternary ground state 
is relevant to the solubility of A (Ref. l37l ). 

TABLE IV: The formation energies AE^ n (Eq. (J27J) for 
Mg-B, B-A, and Mg-A binary ground states (A— Li, Be, Na, 
Mg, K, Ca, Rb, Sr, Cs, Ba). The symbol (+) indicates pos- 
sible ground states never observed experimentally (CaB4 and 
RbB4). The symbol (-) indicates the experimentally observed 
SrgMg38, which was not confirmed to be a ground states by 
first-principle calculations. 





AS at in 

(cV/at.) 


Proto- 
typal 


Space 
groupS 


Pearson 


Ref. 


Li 3 Bi4 


-0.219 


Li 3 Bi4 


I42d(122) 


tI160 


39 


L1B3 


-0.235 


LiB 3 


P4/mbm(127) 


tP20 


40 


Li 8 B 7 


-0.216 


LigBy 


P6m2(187) 


hP15 


32 


Be3Bso 


-0.032 


Be 3 B 50 


P4 2 /nnm(134) 


tP53 


41 


Bei.nB 3 


-0.096 


Bei.nB 3 


P6/mmm(191) 


hPlll 


42 


MgBei 3 


-0.009 


NaZni 3 


Fm3c(226) 


cF112 


43 


NaBis 


-0.059 


NaBi 5 


Imma(74) 


oI64 


44 


Na 3 B 2 o 


-0.070 


Na 3 B 20 


Cmmm(65) 


oC46 


45 


MgB 7 


-0.138 


MgB 7 


Imma(74) 


oI64 


46 


MgB 4 


-0.152 


MgB 4 


Pnma(62) 


oP20 


47 


MgB 2 


-0.151 


A1B 2 


P6/mmm(191) 


hP3 


48 


KB 6 


-0.043 


CaB 6 


Pm3m(221) 


cP7 


49 


CaB 6 


-0.423 


CaB 6 


Pm3m(221) 


cP7 


50 


CaB 4 (+) 


-0.410 


UB 4 


P4/mbm(127) 


tP20 


51 


CaMg 2 


-0.128 


MgZn 2 


P6 3 /mmc(194) 


hP12 


52 


RbB 4 (+) 


-0.163 


UB 4 


P4/mbm(127) 


tP20 




SrB 6 


-0.464 


CaB 6 


Pm3m(221) 


cP7 


53 


Sr 2 Mgi 7 


-0.055 


Th 2 Ni 17 


P6 3 /mmc(194) 


hP38 


54 


Sr 9 Mg3s(-) 


-0.070 


Sr 9 Mg 38 


P6 3 /mmc(194) 


hP94 


55 


Sr 6 Mg 23 


-0.085 


Th 6 Mg 23 


Fm3m(225) 


cF116 


56 


SrM g2 


-0.114 


MgZn 2 


P6 3 /mmc(194) 


hP12 


57 


BaB 6 


-0.421 


CaB 6 


Pm3m(221) 


cP7 


53 


Ba 2 Mgi7 


-0.072 


Zn 17 Th 2 


R3mh(166) 


hR57 


54 


Ba 6 Mg 23 


-0.085 


Th 6 Mg 23 


Fm3m(225) 


cF116 


56 


BaMg 2 


-0.097 


MgZn 2 


P6 3 /mmc(194) 


hP12 


57 



The existence of the ground states for a given binary 
A-B system is based on binary bulk formation energy, 
which, for each phase (j) with stoichiomctry A Xl B XB , is 
determined with respect to pure A and B energies E^t (A) 
and E at (B) as 

A£* t ln (</>) ee £ at (0) - x A Eat(A) - x B E &t (B). (27) 

The results are presented at Table IIVI and in Figure [T] 
For each element A, the reference energy is chosen to be 
the lowest among the pure fee, bee and hep energies^. 
The reference energy for boron is taken to be a-boron 

(Refs. nnmiiiiLil). 

All experimentally observed phases are confirmed ex- 
cept for SrgMg38 (P63/mmc). We also find two possible 




CaMg 2 

1 1 1 1 1 1 1 1 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 



Mg x A A 




1 1 1 1 1 1 1 1 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 

B x A A 



FIG. 1: (Color online) Formation energies A£^ n f° r Mg-B, 
B-A, and Mg-A binary ground states (A=Li, Be, Na, Mg, K, 
Ca, Rb, Sr, Cs, Ba) from Table ITVl The dashed parabola for 
the Li-B system indicates the linear chains phases of Ref. I32I 



new phases, CaB4 and RbB4, both P4/mbm (#127) and 
with UB4 prototype (CaB4 was previously identified in 
ab initio online database^). The formation energies for 
Ba2Mgi7, Ba«Mg23 , BaMg 2 , and CaMg 2 are similar to 
those of Ref. M (reported as -0.079 eV, -0.089 eV, -0.088 
eV, and -0.126 eV, respectively). The small differences 
can be explained on the basis of different GGA pseudopo- 
tcntials (PW91 versus PBE) and different energy cut-offs 
(414 eV versus 360 eV). Our results can not be compared 
with the thermodynamic discussion based on the semi- 
empirical Miedema method reported in Ref. 0, because 
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| MgB 7 / " 


\ U 3 B |4 


Mg-B-Li 


1 MgB 4 'ij-r^- 
1 1 Ai/ 








FIG. 2: (color online) The ternary ground state phase di- 
agrams of Mg-B-A alloys (A=Li, Be, Na). The filled cir- 
cles represent known ground states. The blue stars and red 
triangles represent the substitutional (MgrBseA, Mg3Bi6^4, 
Mg 7 Bi 6 yl) and interstitial (Mg 8 B 56J 4, Mg 4 Bi 6 ,4, Mg 8 Bi 6 ^) 
supercells, respectively. Each supercell has been constructed 
from the most appropriate Mg-B compound (MgB2, MgB4, 
or MgB7). The question mark in the Mg-B-Be system indi- 
cates the existence of new possible ternary ground state(s) 
(Sec. |VI]). 




FIG. 3: (color online) The ternary ground state phase dia- 
grams of Mg-B-A alloys (A=K, Rb, Cs). 



represent the supercells with one interstitial or substitu- 
tional A atom, respectively. In the interstitial case, the 
supercells belong to the lines Mg n B m ^>A, while in the 
substitutional case, the supercells belong to the lines par- 
allel to Mg^4 and intersecting the Mg„B m (constant B 
concentration). 



we include further ground state prototypes other than 
those used in the heat of formations fitting in Ref. [Til . 

The calculated ternary ground state phase diagrams 
for Mg-B-A alloys are depicted in Figs. I2I3I4I Note that 
in each phase diagram, the red triangles and blue stars 



Formation energies numerical results 

The calculated parameters of the model ("raw" /Zj, 7J S 
(Eq. (P)), "true"- e[ i] , e[ s) (Eq. ©) and low-solubility 
-^soi' ^"soi (Eqs. ([3])) formation energies) are presented 
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FIG. 4: (color online) The ternary ground state phase dia- 
grams of Mg-B-A alloys (A=Ca, Sr, Ba) . 



in Table E] and in Fig. G3 

The values of the low-solubility energies reported in 
Table |V] suggest the following, (a) Only E^J-^ for sub- 
stitutional Be in MgB7 was found to be negative, indi- 
cating the existence of unknown ground state(s) within 
the triangle MgBy ^603650 <->Bei.nB3. This fact is 
summarized by the question mark in Fig. [2j (b) Gen- 
erally, the substitutional systems have lower formation 
energy than the interstitial systems. Exceptions are Be, 
Cs, and Ba in MgB2; (c) Although oscillating, the forma- 
tion energies tends to increase with the element number, 
as shown in Fig. Regular oscillations are observed 



TABLE V: Interstitial (i) and substitutional (s) quantities for 
alkali and alkaline earth solute elements A in MgB2, MgBzt, 
and MgB7 compounds.: "raw" formation energies m, 71 s (Eq. 
CO), "true" formation energies E ( A s) (Eq. ©) and low- 

solubility formation energies E^-^, E^?\ (Eqs. ©). In each 
case, the presented value is the lowest one among all the pos- 
sible interstitial or substitutional positions. 







Mg 


B 2 


Mg 


B 4 


Mg 


B 7 


A 


(eV) 


(i) 


(s) 


(i) 


(s) 


(i) 


(s) 


Li 


Mi, Ms 


0.542 


-0.111 


0.234 


-0.109 


-0.959 


-0.055 




Z7>(i>s) 


2.444 


0.310 


2.136 


0.312 


0.943 


0.366 




E (x,*) 
sol 


2.716 


0.574 


2.654 


0.686 


1.657 


0.506 


Be 




-2.665 


-0.876 


-0.506 


-0.710 


-1.577 


-2.100 




E u 


1.042 


1.350 


3.201 


1.516 


2.130 


0.126 






1.054 


1.21 


3.201 


1.256 


2.181 


-0.557 


Na 


Mi, Ms 


5.623 


1.839 


2.409 


1.285 


1.999 


0.869 




Ei(i> s ) 


6.932 


1.666 


3.718 


1.113 


3.308 


0.697 






6.944 


1.526 


3.718 


0.802 


3.308 


0.147 


K 




7.022 


4.674 


2.984 


3.001 


3.453 


2.077 




E u 


8.061 


4.232 


4.023 


2.559 


4.492 


1.635 






8.073 


4.092 


4.023 


2.248 


4.492 


0.838 


Ca 




2.833 


-0.295 


1.042 


-1.107 


1.251 


-1.941 




Ei(i>") 
\.\ 


4.751 


-0.732 


2.960 


-0.671 


3.169 


-1.504 




r-i(l.s) 


6.364 


1.747 


4.993 


1.218 


5.453 


0.358 


Rb 


Un . U 


9.463 


6.572 


3.733 


3.957 


4.751 


3.312 






10.393 


6.021 


4.663 


3.406 


5.681 


2.761 






10.406 


5.981 


4.862 


3.460 


6.047 


2.477 


Sr 


Mi i M s 


5.514 


1.761 


1.334 


0.221 


2.149 


-0.960 




f ^ 


7.138 


1.618 


2.958 


0.364 


3.773 


-0.817 






9.042 


3.802 


5.282 


2.544 


6.348 


1.331 


Cs 




4.663 


8.436 


4.559 


5.064 


6.488 


5.131 




E M 


5.523 


7.815 


5.419 


4.444 


7.349 


4.510 






5.536 


7.675 


5.419 


4.134 


7.349 


3.411 


Ba 




2.450 


3.959 


1.365 


1.084 


3.118 


0.335 




E M 


4.373 


3.518 


3.288 


1.525 


5.041 


0.777 






5.971 


5.99 


5.306 


3.400 


7.309 


2.624 



for all substitutionals systams starting with Na. The 
lower and higher boundaries of such correspond to alka- 
line earth and alkali elements, respectively, (d) Na, Ca, 
and Li substitutionals in MgBj have the lowest formation 
energies (0.147 eV, 0.358 eV, and 0.506 eV, respectively), 
resulting in high solubility (Sec. IVIIj) . (e) In general for 
substitutionals systems, the higher boron contents the 
lower formation energy (except Li and Be), (f) The cor- 
responding "true" E A ' S ^ and low-solubility E^^(A) for- 
mation energies generally demonstrate a similar behavior 

(s) 

as functions of A. Exceptions are E x k ' < for A=Ca, and 

Sr and E^^A) < for A=Bc. The difference is due to 
the consideration of all the known ground states for the 
determination of E^J and not only pure Mg, B and A 
(as for E^ s) ). 
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A- Li Be Na K Ca Rb Sr Cs Ba 



FIG. 5: (Color online) Representation of the formations en- 
ergies from Table |V] Dashed (green) and dotted (red) lines 
connect the lowest formation energies for alkali and alkaline 
earth elements, respectively. 



For Li-Mg-B, in Ref. [I?] the formation energies were 
obtained from first principles for interstitial and substi- 
tutional Li in boron-rich MgB2 (3.00 eV and 0.27 eV, 
respectively). The results of Ref. [ij] are different form 
our low-solubility formation energies -Ego 1 (Li)=2.72 eV 

and _E^i( Li )= - 57 cV ( scc Table Ei because only the 
hexagonal a-LIB phase was considered as a ground state 
in Ref. [13 and the other boron-rich phases were reported 
to be unstable (the a-LiB was refined in Ref. [3^ by vari- 
ational minimization of the Li and B concentration as 
LigBy). However, in our high-throughput framework, we 
found L1B3 and Li3Bi4 to be stable^, and the appro- 
riate analysis of the formation energies should be done 
with respect to the triangle MgB2^->LiB3^Mg as shown 
in Fig. [2] In the low-solubility limit of Li, since the su- 
percell concentration is close to MgB2, the lever rule of 
phase decomposition of the mixtures does not cause sig- 
nificant errors if the energies are calculated with respect 
to other references, and therefore the results of Ref. [13 
are somehow similar. 

For Na-Mg-B, it can be shown that the substitutional 
Na formation energy in magnesium-rich MgB2 (1.74 eV) 
obtained in Ref. [IT] is equivalent to our "true" formation 



energy for substitutional Na in MgB 2 (E£> = 1.67 cV, 
see Table |Vj . The observed small numerical difference 
between two results can be attributed to the different 
energy cutoffs for the basis set (414 eV versus 312 eV)££. 

VII. SOLUBILITY RESULTS 

Solubility results are presented in Table I VII (the low- 
solubility approximation values are shown in brackets). 
We report only values larger than 10 -6 (Li, Na, Ca, and 
K). The values on Table IVT1 suggest the following, (a) 
Na is the only element which has a substantial solubil- 
ity: Na in MgB 7 is - 0.5-1 % at T = 650 - 1000 K. (b) 
The substitutional solubility of Be in MgBy can not be 
determined because the corresponding formation energy 

^soi ( see Eq- 0) is negative implying that the ground 
state list is not complete. Thus the approach of Sec. 
fVl is not applicable, (c) In agreement with their lowest 

formation energies 

£^1 (Table Eland Fig. [5]), substitu- 
tional Li, Na, and Ca experience the largest solubilities 
among all studied systems; (d) Due to their comsiderably 
high formation energies (Table El and Fig. [5J), all 

the investigated interstitial systems have negligible sol- 
ubilities, (e) The low-solubility approximation and the 
general theory agree within 5%. 

The calculated negligible solubilities of Be, Na, 
Ca in MgB2 agree with previous experimental 
observations^ ^ 11,12^3 ^ e n0 ^ nnc j an y g jg n 

of nonvanishing solubilities of Rb, Cs and Ba in the 
MgB2 alloy even at high annealing temperature in 
disagreement with the reported experimental data of 
Ref. 1. However, our negligible Rb and Cs solubility 
in bulk MgB 2 are in agreement with the experimental 
conclusion made in Ref. [13 where the author suggest 
that Rb and Cs dopants most likely segregate in grain 
boundaries. 

The obtained low solubilities of Li and Sr in MgB2 
differ from experimental vales (solubility up to 30% in 
Refs. iillilllll). This discrepancy can be attributed 
to segregation of Li and Sr in the grain boundaries as 
was concluded for Rb and Cs in Ref. [l(| In particular, 
the value of E^'^3.8 cV for Sr is too large for non- 
negligibllc bulk solubility. For Li in MgB2 the calculated 

(s) 

formation energy i?,: ol =0.57 eV is not very large and it 
is the smallest among the investigated solutes in MgB2- 
In Ref. [13, qualitative conclusions about high solubility 
of Li and low solubility of Na in MgB2 were made based 
on a large difference between the corresponding forma- 
tion energies. Our numerical results demonstrate that, 
despite much lower formation energy of Li in comparison 
with Na in MgB2, the Li-solubility is still very low (Sec. 

ED- 

It should be emphasized that our model is in thermo- 
dynamical equilibrium which can be difficult to reach at 
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TABLE VI: Calculated equilibrium interstitial x^ (i) and 

substitutional %j^ s ) ( s ) A-solute concentrations at tempera- 
tures T =300, 650, and 1000 K, in MgB 2 , MgB 4 , and MgB 7 
compounds. The values obtained within the low-solubility 
approximation are reported in brackets. Zeros are used for 
values smaller than 10 -6 . For Be "AE < 0" indicates that 
the corresponding formation energy EQ, (see Eq. ©) is neg- 
ative implying that the ground state(s) list is not complete 
and the approach of Sec. [V]is not applicable. 



X 


(i) 


MgB 2 


MgB 4 


MgB 7 


A 


T (K) 


i 


s 


i 


s 


i 


s 


Li 


300 
























650 





1.3x10" 


-5 











9.9xl0 -6 








(1.2x10" 


■ 5 ) 








(9.9xl0~ 6 ) 




1000 





4.5x10" 


-4 





6.9xl0" 5 





2.6xl0"~ 4 








(4.3x10" 


- 4 ) 




(7.0xl0 -5 ) 




(2.6xl0 -4 ) 


Be 


300 



















AE < 




650 



















AE < 




1000 



















AE < 


Na 


300 



















2.1xl0~ 4 


















(2.1xl0 -4 ) 




650 



















4.4xl0" 3 


















(4.2xl0 -3 ) 




1000 













1.8xl0~ 5 





l.OxlO" 2 














(1.8xl0 -5 ) 




(0.96xl0 -2 ) 


K 


ouu 
























650 
























1UUU 



















3.6xl0~ 6 




OUU 
























650 



















l.lxl0~ 4 


















(l.OxlO -4 ) 




1000 



















l.OxlO" 3 


















(0.96xl0 -3 ) 


Rb 


300 
























650 
























1000 






















Sr 


300 
























650 
























1000 






















Cs 


300 
























650 
























1000 






















Ba 


300 
























650 
























1000 























low temperatures and the experimental equilibrium solu- 
bility tends usually to be overestimated. In fact, the for- 
mation of metastable and/or unstable states which are 
subsequently frozen at low temperatures, can make solu- 
bility measurements very challenging. In such scenarios, 
the measured solubility may correspond to spinodal con- 
centration rather then actual binodal concentration or 
simply characterize the frozen out of equilibrium solubil- 
ity remaining from the initial specimen preparation at 
higher temperature. Besides, the segregation of defects 



into grain-boundaries, especially in multicrystallinc sam- 
ples prepared through non optimal cooling dramatically 
affect the amount of frozen defects and solutes. 

Although in our study we did not perform an extensive 
search over the configurational space of ternary alloy, if 
a new ternary phases were present near MgB2, the so- 
lute atoms would concentrate and nucleate such phases, 
and this may be misinterpreted as a high solubility in 
MgB2 phase. Furthermore, a new equilibrium ternary 
phase would result in an increase of impurities formation 
energies and, correspondingly, in a decrease of solubility. 
Numerical approximations in the first-principles arc not 
expected to affect the values of solubility: i.e. tipically 
for SE so i ~ 30 mcV/atom, Sx^ ~2x 10~ 4 . Lattice vi- 
brations and solute-solute interactions are neglected be- 
cause they are not expected to play important roles at 
low temperatures and low solute concentrations. 



VIII. CONCLUSIONS 

In the present paper, we present an approach to study 
the solubilities in ternary alloys. The advantage of the 
approach is in taking into account all known ternary 
ground states rather than just pure solids. Based on 
the approach, we propose an analytical low-solubility ap- 
proximation that can be used for high-throughput calcu- 
lations of solubilities in alloys. 

Combining the developed approach with first principle 
calculations, we have determined the formation energies 
and solubilities of alkali (Li, Na, K, Rb, Cs) and alkaline 
earth (Be, Ca, Sr, Ba) metals in the Mg-B system. It is 
found that the considered metals have low solubilities in 
the boron-rich Mg-B alloy. Substitutional Na, Ca, and Li 
experience the the largest solubilities, with Na in MgB 7 
reaching 0.5-1% at T = 650 - 1000 K. All the considered 
interstitial scenarios leed to negligible solubilities. The 
solubility of Be in MgB 7 can not be determined with our 
model because the corresponding low-solubility forma- 
tion energy is negative implying that the existing ground 
states list must be augmented through a more extensive 
search over the configurational space. 

We also present a high-throughput search of ground 
states in binary Mg-B, Mg-^4, and B-A alloys (A=Li, 
Be, Na, Mg, K, Ca, Rb, Sr, Cs, Ba). Ternary phase 
diagrams Mg-B- A are constructed based on of the de- 
termined phases. SrgMg3g is not an equilibrium ground 
state despite its high temperature validations. Two new 
ground states CaB 4 and RbB 4 are found. 
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